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ABSTRACT 

We present the sub-millimeter spectra from 450 GHz to 1550 GHz of eleven nearby active galaxies 
observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) onboard Herschel. Wc 
detect CO transitions from J up = 4 to 12, as well as the two [C I] fine structure lines at 492 and 
809 GHz and the [Nil] 1461 GHz line. We used radiative transfer models to analyze the observed 
CO spectral line energy distributions (SLEDs). The FTS CO data were complemented with ground- 
based observations of the low-J CO lines. We found that the warm molecular gas traced by the 
mid-J CO transitions has similar physical conditions (tih 2 <~10 3 ' 2 — 10 3 ' 9 cm -3 and T^ n ^300-800 K) 
in most of our galaxies. Furthermore, we found that this warm gas is likely producing the mid-IR 
rotational H2 emission. We could not determine the specific heating mechanism of the warm gas, 
however it is possibly related to the star- formation activity in these galaxies. Our modeling of the 
[C I] emission suggests that it is produced in cold (Xkin < 30 K) and dense (nn 2 > 10 3 cm -3 ) molecular 
gas. Transitions of other molecules arc often detected in our SPIRE/FTS spectra. The HF J = 1 — 
transition at 1232 GHz is detected in absorption in UGC 05101 and in emission in NGC 7130. In the 
latter, near-infrared pumping, chemical pumping, or collisional excitation with electrons are plausible 
excitation mechanisms likely related to the AGN of this galaxy. In some galaxies few H2O emission 
lines are present. Additionally, three OH + lines at 909, 971, and 1033 GHz are identified in NGC 7130. 
Keywords: galaxies: active - galaxies: ISM - galaxies: nuclei - galaxies: Seyfcrt 



1. INTRODUCTION 

The sub-millimeter (sub-mm) is one of the few as- 
tronomical spectral regions that has not been fully 
explored so far. In part, this is because ground- 
based sub-mm observations arc difficult due to the 
low transparency of the atmosphere. In consequence, 
only a small number of bright sources have been ob- 
served in a few atmospheric windows (350, 450, and 
850 /Ltm). The SPIRE Fourier transform spectrome- 
! ter (SPIRE/FTS; Griffin et al. 2010; Naylor et al. 2010; 
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Swinyard et al. 2010) onboard the Herschel Space Ob- 
servatory (Pilbratt et al. 2010) fills most of this gap cov- 
ering the spectral range between 450 and 1440 GHz (210 
and 670 /im) and provides an important benchmark for 
high-redshift studies with the Atacama Large Millime- 
ter/submillimeter Array (ALMA). 

Besides the continuum arising from cold dust, 
in the sub-mm range there are several molecular 
and atomic spectral lines that probe the different 
phases of the interstellar medium (ISM) in galax- 
ies (e.g., van der Werf et al. 2010; Panuzzo et al. 2010; 
Rangwala et al. 2011). In the SPIRE/FTS spectral 
range the most prominent spectral feature is the CO lad- 
der from J up = 4 to 13. CO is a good tracer of the molec- 
ular gas since it is one of the most abundant molecules 
in the ISM. In particular, these mid-J CO lines originate 
in warm molecular gas (their upper-level energies range 
from 55 to 500 K above the ground state) that can be ex- 
cited by ultraviolet (UV) photons in photon-dominated 
regions (PDRs; e.g., Wolfire et al. 2010), X-rays in X-ray 
dominated regions (XDRs; e.g., Meijerink et al. 2006) 
or by shocks (e.g., Flower & Pineau Des Forets 2010). 
The analysis of the CO spectral line energy distri- 
bution (SLED), including low- and mid-J CO tran- 
sitions, is crucial for understanding the conditions in 
the cold and warm molecular gas and the dominant 
heating mechanism in galaxies (Rangwala et al. 2011; 
Kamenetzky et al. 2012; Spinoglio et al. 2012). 

In addition to the CO ladder, two [C I] fine struc- 
ture lines are present in the sub-mm range. PDR models 
predict that these lines originate in the transition re- 
gion between C + and CO (Tielens & Hollenbach 1985; 
Kaufman et al. 1999). However, several mechanisms, 
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Table 1 

Sample of Local Active Galaxies 



Galaxy 


R.A. 


Decl. 


cz a 


Dist. b 


Spect. Class c 


Rcf. d 


log I/IR 6 




Ref.s 




(J2000.0) 


(J2000.0) 


(km s" 1 ) 


(Mpc) 






(£©) 


(10 41 ergs" 1 ) 




NGC 1056 


02 42 48.3 


+28 34 27 


1545 


32.2 


H II 


1 


10.3 






UGC 05101 


09 35 51.6 


+61 21 11 


11802 


163 


LINER 


1 


12.0 


75 


7 


NGC 3227 


10 23 30.6 


+19 51 54 


1157 


14.4 


Syl.5 


2 


9.9 


1.9 


8 


NGC 3982 


11 56 28.1 


+55 07 31 


1109 


20.6 


Syl-9 


2 


10.0 


<0.025 


9 


NGC 4051 


12 03 09.6 


+44 31 53 


700 


14.6 


Syl.2 


2 


10.0 


1.3 


10 


NGC 4151 


12 10 32.6 


+39 24 21 


995 


12.1 


Syl.5 


2 


9.8 


8.2 


10 


NGC 4388 


12 25 46.7 


+12 39 44 


2524 


18.2 


Syl.9/HBLR 


2, 3 


10.2 


3.0 


11 


IC 3639 


12 40 52.8 


-36 45 21 


3275 


42.5 


Sy2/HBLR 


1, 4 


10.8 


0.17 


12 


NGC 7130 


21 48 19.5 


-34 57 04 


4842 


66.0 


Sy2 


1 


11.4 


0.83 


13 


NGC 7172 


22 02 01.9 


-31 52 11 


2603 


33.9 


Sy2 


5 


10.4 


29 


10 


NGC 7582 


23 18 23.5 


-42 22 14 


1575 


20.6 


H Il/Syli 


1, 6 


10.8 


2.0 


14 



References. — (1) Yuan et al. (2010); (2) Ho et al. (1997); (3) Young et al. (1996); (4) Heisler et al. (1997); 
(5) Lumsden et al. (2001); (6) Reunanen et al. (2003); (7) Imanishi et al. (2003); (8) Gondoin et al. (2003); 
(9) Guainazzi et al. (2005); (10) Brightman & Nandra (2011); (11) Cappi et al. (2006); (12) Guainazzi et al. 
(2005); (13) Levenson et al. (2005); (14) Piconcelli et al. (2007). 



Note. — (a) Heliocentric velocity from NASA/IPAC extragalactic database (NED). (i>) Distance from NED. 
^ Nuclear activity classification. Syli indicates that broad emission lines are observed in the near-IR. HBLR 
indicates the detection of a hidden broad line region. ^ Reference for the nuclear activity classification. ' e ' 
8-1000 /im infrared luminosity from Sanders et al. (2003) scaled to our adopted distance. ^' X-ray 2-10 keV 
luminosity. ^ Reference for the X-ray luminosity. 



such as turbulent diffusion, non chemical equilibrium, or 
intense X-ray radiation, can influence the C° distribution 
and abundance (Papadopoulos et al. 2004 and references 
therein) . 

To date, there are relatively few published SPIRE/FTS 
spectra of extragalactic objects (van der Werf et al. 
2010; Panuzzo et al. 2010; Rangwala et al. 2011; 
Kamcnctzky et al. 2012; Spinoglio et al. 2012; 
Meijerink et al. 2013); in these, however, molecules 
such as H2O, OH + , and HF arc often detected. In 
addition, other molecules like HCN, CH + , and H2 + 
have also been detected. 

In this paper we present for the first time the com- 
plete sub-mm spectra of a sample of eleven active 
galaxies (Section 2). The sub-mm spectra are comple- 
mented by ground-based IRAM 30 m observations of the 
CO J = 1 — and J = 2 — 1 transitions in three galax- 
ies. The SPIRE and IRAM 30 m observations and data 
reduction are described in Section 3. We use the es- 
cape probability approximation to model the CO and C° 
emission and study how they trace the warm and cold 
molecular gas in these galaxies (Sections 4 and 5). In 
two galaxies we detect the HF J = 1 — transition, one 
in emission and one in absorption. We discuss several 
possible origins for this line in Section 6. Finally, a brief 
discussion about the H2O and OH + lines detected in a 
few galaxies is presented in Section 7. The [N II] line 
at 205 /im (1461GHz) is one of the brightest far-IR lines 
in galaxies (e.g., Wright et al. 1991), however it has not 
been systematically explored by previous observatories. 
In a forthcoming paper we will study the emission of the 
[N II] 205 /im line and other far-IR atomic fine structure 
emission lines in our sample of galaxies. Throughout this 
paper we assumed Hq = 73kms~ 1 Mpc -1 . 

2. SAMPLE OF LOCAL ACTIVE GALAXIES 

We study a sample of eleven local (d = 12 — 160 Mpc) 
active galaxies drawn from the 12 /im Galaxy Sam- 



ple (Rush et al. 1993). They were selected because of 
their high IRAS 100 /im flux and to include the differ- 
ent nuclear activity classes in similar proportions (three 
Seyferts type 1.2 and 1.5, two type 1.9 Seyfert, three type 
2 Seyfert, one low-ionization nuclear emission line region 
[LINER], and two H II galaxies). Their total infrared 
(IR) luminosities (Lir) range from 10 9 8 to 1O 12 O L , 
and the median value is 10 10 ' 3 Lq. In our sample, only 
UGC 05101 and NGC 7130 are classified as luminous 
IR galaxies (LIRGs; Lm > lO n L ), so in these galax- 
ies, besides the AGN activity, intense star-formation is 
expected (Petric et al. 2011; Alonso-Herrero et al. 2012). 
The galaxies in our sample are presented in Table 1. 

3. OBSERVATIONS AND DATA REDUCTION 

All the galaxies in our sample were observed with 
the SPIRE/FTS and the SPIRE photometer onboard 
Herschel. Additionally, we observed the CO J = 1 — 
and J = 2 — 1 rotational transitions in UGC 05101, 
NGC 3227, and NGC 3982 with the IRAM 30 m tele- 
scope. The SPIRE photometric data will be analyzed 
in detail in a forthcoming paper, but the data reduction 
is summarized here because the images are used to cal- 
culate the source-beam coupling factor and the scaling 
factor for the flux calibration of the SPIRE/FTS spec- 
tra. 

3.1. SPIRE/FTS Spectroscopy 

We collected Herschel SPIRE/FTS spectroscopic ob- 
servations of eleven nearby active galaxies through a 
guaranteed time project (PI: L. Spinoglio) and a guaran- 
teed time key project (PI: C. D. Wilson). The integra- 
tion times of the observations were ~5100s, achieving a 
sensitivity of 0.2-0.3 Jy in the continuum and a 3c de- 
tection limit of ~1 x 10~ 14 erg cm" 2 s" 1 for an unresolved 
spectral line. 

The SPIRE/FTS has two bolometer arrays: the spec- 
trometer short wavelength (SSW; 194-324 /mi, 925- 
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Table 2 

Log of Herschel Observations 



Galaxy 


SPIRE/FTS Spectrometer 
Observation ID 


SPIRE Photometer 
Observation ID 


NGC 1056 


1342204024 


1342226630 


UGC 05101 


1342209278 


1342204962 


NGC 3227 


1342209281 


1342197318 


NGC 3982 


1342209277 


1342186862 


NGC 4051 


1342209276 


1342210502 


NGC 4151 


1342209852 


1342188588 


NGC 4388 


1342210849 


1342211416 


IC 3639 


1342213381 


1342202200 


NGC 7130 


1342219565 


1342210527 


NGC 7172 


1342219549 


1342209301 


NGC 7582 


1342209280 


1342210529 



1545 GHz) and the spectrometer long wavelength (SLW; 
316-672 lira, 446-948 GHz). The beam full-width half- 
maximum (FWHM) of the SSW bolometers is ~18", ap- 
proximately constant with frequency. On the contrary, 
the beam FWHM of the SLW bolometers varies between 
~30" and ~42" with a complicated dependence on fre- 
quency (Swinyard et al. 2010). They sample sparsely a 
field of view (FOV) of ~2', with the central bolometer 
of the arrays centered at the nuclei of the galaxies (see 
Figure 1). At the distances of our galaxies 30" corre- 
sponds to physical scales of 2 to 20kpc. All the galax- 
ies were observed at high spectral resolution (unresolved 
line FWHM 1.45 GHz, u/Av ~ 300- 1000) in the single 
pointing mode. The observation IDs arc listed in Table 
2. 

We processed the raw data with the Herschel inter- 
active pipeline environment software (HIPE) version 9 
(Ott 2010). We used the standard reduction scripts 
(Fulton ct al. 2010). These scripts process the raw 
SPIRE/FTS timelines and produce a flux calibrated 
spectrum for each bolometer. Briefly, first the bolometer 
timelines are corrected for instrumental effects and cos- 
mic ray hits, and the interferograms are created. The 
baselines are removed and the phase of the interfero- 
grams corrected before a Fourier transform is applied to 
obtain the spectra. These spectra are dominated by the 
telescope (and instrument) thermal emission that must 
be removed to obtain the source spectrum. 

The uncertainty in the telescope emission model used 
by HIPE is 1-2 Jy in the continuum level (SPIRE Ob- 
server's Manual 12 ), which is comparable to the fluxes of 
our galaxies. Therefore, the residual telescope emission 
has to be subtracted before applying the point source cor- 
rection to the spectra. The optical angular sizes of our 
galaxies are several arcmin, but the nuclear far-IR emis- 
sion appears compact at the angular resolution of SPIRE 
(only UGC 05101 and NGC 7130 are almost point-like 
at the SPIRE/FTS angular resolution, see Figure 1). 
Moreover, far-IR emission is detected only in the cen- 
tral bolometer of the SPIRE/FTS arrays (FWHM-18" 
and 30-40" for the SSW and SLW arrays, respectively). 
Thus, we estimated the residual telescope emission av- 
eraging the spectra of the off-axis bolometers. For some 
galaxies, some of the off-axis detectors contain extended 
galactic emission (weak continuum and [Nil] 1461 GHz 

12 Available at http://herschel.esac.esa.int/Docs/SPIRE/ 
html/spire_om. html 



emission). We excluded those detectors before averag- 
ing. 

Next we applied the point source correction 13 to the 
central bolometer spectra to obtain the flux calibrated 
spectra of the galactic nuclei. To correct for the large 
beam size variations in the SLW spectra we calculated 
the source-beam coupling factor, 77. We define 77 as the 
ratio between the flux measured with a beam of size 9 
and that measured with a 30" beam (the minimum SLW 
beam size). Similar to Kamenetzky et al. (2012) and 
Panuzzo et al. (2010), we convolved the 250 /mi SPIRE 
maps using Gaussian profiles to obtain the coupling fac- 
tor as a function of the beam size, r/(9). Then with the 

1 

SLW beam size dependence on the frequency , 0(f), we 
obtained r](y) and we divided the SLW spectra by this 
function. The values of rj(v) for these galaxies vary be- 
tween 1.0 and 1.4. By applying this correction we assume 
that the spatial distributions of the emission lines in the 
SPIRE/FTS spectra are comparable to that of the cold 
dust traced by the 250 /im emission. 

Finally, we scaled the SSW and SLW spectra to match 
the SPIRE photometry (see Section 3.2). For the SSW 
we scaled the spectra to match the 250 /mi (1200 GHz) 
flux measured in a 18" diameter aperture. Similarly, we 
scaled the SLW spectra to match the flux measured in the 
350 /im (856 GHz) image in a 30" diameter aperture. We 
used the 350 /tm SPIRE image because for most of our 
sources the continuum at 500 /im (600 GHz) is below the 
detection limit of the FTS data. Point source aperture 
corrections are applied to the SPIRE photometry (see 
Section 3.2). The scale factors vary between 0.85 and 
1.30 for the SSW spectra and between 0.95 and 1.20 for 
the SLW spectra. For NGC 1056, NGC 4051, NGC 4151, 
IC 3639, and NGC 7172 the 350 /tm continuum is not 
detected in the FTS spectra. Thus for these galaxies 
the uncertainty in the absolute flux of the emission lines 
in the SLW range is about 20 % larger. The final flux 
calibrated spectra of the galaxies are shown in Figure 2. 

We measured the line intensities by fitting a sine pro- 
file 14 to the emission lines. To obtain the local continuum 
we used a linear fit in a ^10 GHz range around the line 
position. When two lines were close in frequency (e.g., 
CO J = 7 - 6 at 807 GHz and [C I] 809 GHz) the two lines 
were fitted simultaneously. In the SPIRE/FTS spectra 
of our galaxies we find CO lines from J up = 4 to J up = 
12, the [Cl]492GHz, [Cl]809GHz, and [Nil] 1461 GHz 
atomic lines, and some water lines. The HF J = 1 —0 ro- 
tational transition at 1232 GHz is detected in two sources 
and in one source we detect three OH + transitions. The 
fluxes and la uncertainties of the detected lines and the 
3<7 upper limits for the undetected lines are given in Ta- 
bles 3 and 4. 

3.2. SPIRE Photometry 

We obtained Herschel SPIRE imaging observations for 
seven of our galaxies through a guaranteed time project 
(PI: L. Spinoglio). The rest of the galaxies in our sam- 
ple were observed as part of the Very Nearby Galaxies 

13 Taken from the HIPE version 9 calibration product 
BeamParam_HR_unapod_nominal_20050222. It provides the point 
source correction and the beam size as a function of frequency 
with a 0.3 GHz sampling. 

14 The sine function represents well the SPIRE/FTS instrumen- 
tal line shape (see the SPIRE Observer's Manual). 
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Table 3 

SPIRE/FTS Line Fluxes 



Tr&iisitioii 








Fluxes (10~ 15 


erg cm ^ s ^ ) 








(GHz) 


NGC 1056 


UGC 05101 


NGC 3227 


NGC 3982 


NGC 4051 


NGC 4151 


CO J = 4-3 


461.041 


<15.8 




29.0 ± 3.7 


18.1 ± 2.7 


<13.0 


<9.5 


12 CO J = 5-4 


576.268 


10.1 ± 2.0 


15.6 ± 3.8 


33.4 ± 3.9 


12.9 ± 2.2 


<11.2 


<8.0 


12 CO J = 6-5 


691.473 


6.0 ± 1.5 


16.9 ± 2.3 


25.2 ± 2.9 


7.3 ± 1.6 


9.6 ± 1.6 


<4.1 


12 CO J = 7-6 


806.652 


6.1 ± 1.4 


11.1 ± 1.8 


18.1 ± 2.3 


4.7 ± 1.3 


7.3 ± 1.4 


3.6 ± 1.1 


12 CO J = 8-7 


921.800 


<8.0 


11.6 ± 2.6 


23.7 ± 2.9 


<10.1 


7.6 ± 1.6 


5.3 ± 1.3 


12 CO J = 9-8 


1036.912 


<13.6 


<11.7 


14.6 ± 2.3 


<12.1 


<15.1 


<8.1 


12 CO J = 10-9 


1151.985 


<13.3 


13.1 ± 2.5 


15.7 ± 2.4 


<12.9 


<12.2 


<9.8 


12 CO J = 11-10 


1267.015 


<11.6 


9.7 ± 2.9 


13.9 ± 2.4 


<13.6 


<17.8 


<9.0 


12 CO J = 12-11 


1381.995 


<12.1 


<8.7 


<14.2 


<13.0 


<15.0 


<9.0 


12 CO J = 13-12 


1496.923 


<16.7 


<11.5 


<14.4 


<15.6 


<19.7 


<11.4 


p-H 2 2 n -2 2 


752.033 


<5.8 


8.5 ± 1.8 


<7.5 


<5.3 


<5.3 


<3.1 


p-H 2 2 02 -lii 


987.927 


<18.2 


<11.0 


<14.4 


<15.3 


<20.7 


<10.9 


O-H2O 3l2~3o3 


1097.365 


<11.6 


<9.6 


<12.9 


<11.0 


<12.2 


<8.1 


P-H2O ln-Ooo 


1113.343 


<13.0 


<11.4 


<13.7 


<10.6 


<12.2 


<8.4 


0-H 2 321~3l2 


1162.912 


<14.1 


10.8 ± 2.4 


<11.2 


<12.0 


<12.3 


<10.2 


p-H 2 2 20 -2n 


1228.789 


<11.9 


<12.9 


12.0 ± 3.2 


<11.2 


<16.4 


<9.5 


[C I] 3 Pi- 3 P 


492.161 


<12.2 


<16.0 


22.7 ± 4.3 


<13.2 


<11.3 


<6.9 


[C I] 3 P 2 - 3 Pi 


809.342 


7.7 ± 1.5 


15.8 ± 2.1 


46.6 ± 1.9 


10.0 ± 1.5 


7.8 ± 1.3 


8.4 ± 1.0 


[N II] 3 Pi- 3 P 


1461.128 


50.5 ± 3.4 


37.9 ± 4.4 


38.6 ± 4.3 


88.1 ± 3.7 


25.0 ± 4.3 


21.2 ± 3.1 



Note. — Measured line fluxes and la statistical uncertainties. For the non-detections we state the 
3a upper limits. ^ Due to the higher redshift of UGC 05101 the CO J = 4 — 3 transition is outside the 
SPIRE/FTS range. 



Table 3 

SPIRE/FTS Line Fluxes - Continued 



Transition 


t'rest 




Fluxes (10 15 ergcm 


- 2 s^) 






(GHz) 


NGC 4388 


IC 3639 


NGC 7130 


NGC 7172 


NGC 7582 


12 CO J = 4-3 


461.041 


23.9 ± 2.9 


12.1 ± 1.9 


35.3 ± 3.9 


<53.9 


69.3 ± 7.2 


12 CO J = 5-4 


576.268 


19.1 ± 2.4 


<10.3 


36.6 ± 3.9 


<43.7 


78.0 ± 8.0 


12 CO J = 6-5 


691.473 


15.4 ± 1.9 


<5.4 


27.7 ± 3.0 


21.2 ± 2.8 


82.3 ± 8.3 


12 CO J = 7-6 


806.652 


13.1 ± 1.7 


<4.8 


27.8 ± 2.9 


<19.0 


65.6 ± 6.7 


12 CO J = 8-7 


921.800 


14.1 ± 1.8 


<7.3 


22.2 ± 2.6 


<27.9 


63.1 ± 6.5 


12 CO J = 9-8 


1036.912 


9.7 ± 1.6 


<9.6 


25.2 ± 2.8 


<9.3 


54.6 ± 5.7 


12 CO J = 10-9 


1151.985 


<8.1 


<9.2 


18.8 ± 2.4 


<10.4 


35.1 ± 3.8 


12 CO J = 11-10 


1267.015 


<10.8 


<11.2 


12.9 ± 1.8 


<12.3 


26.4 ± 3.0 


12 CO J = 12-11 


1381.995 


<8.0 


<8.6 


10.3 ± 1.6 


<9.3 


10.9 ± 1.9 


12 CO J = 13-12 


1496.923 


<12.2 


<11.4 


<12.4 


<11.9 


<22.4 


p-H 2 2 n -2 2 


752.033 


<5.8 


<4.6 


6.6 ± 0.9 


<18.1 


<9.2 


p-H 2 2 02 -lii 


987.927 


<11.0 


<11.7 


18.6 ± 2.0 


<12.1 


<16.9 


O-H2O 3l2-3o3 


1097.365 


<8.1 


<9.6 


12.2 ± 1.5 


<9.0 


<14.0 


p-H 2 ln-0,,0 


1113.343 


<8.3 


<8.8 


<9.5 


<9.6 


<11.8 


O-H 2 321-3l2 


1162.912 


<16.7 


<8.2 


<12.8 


<18.7 


13.9 ± 2.1 


P-H2O 2 20 -2ii 


1228.789 


<11.0 


<10.0 


12.8 ± 1.1 


<9.9 


<13.1 


[C I] 3 Pi- 3 P 


492.161 


14.5 ± 2.0 


<11.7 


17.9 ± 2.2 


38.2 ± 5.2 


34.4 ± 2.8 


[C I] 3 P 2 - 3 Pi 


809.342 


26.7 ± 0.7 


5.8 ± 0.7 


28.2 ± 1.0 


45.0 ± 4.0 


75.9 ± 2.0 


[N II] 3 Pi- 3 P 


1461.128 


47.4 ± 1.9 


45.5 ± 2.1 


123.0 ± 1.8 


120.1 ± 4.1 


191.8 ± 5.4 
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Figure 1. Footprint of SPIRE/FTS FOV plotted over the Herschel /SPIRE 250 /im images of our sample of Seyfert galaxies. The solid 
and dashed circles mark the positions of the SSW and SLW bolometers, respectively. Only the functioning detectors in the unvignetted 
FOV are plotted. North is up and East to the left. The images are shown in a linear gray scale. The black line represents 3kpc at the 
distance of the galaxy. 
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Figure 2. Observed SPIRE/FTS spectra of our sample. The black and gray lines are the SSW and SLW spectra, respectively. Notice the 
overlap between the two spectra in the 960 and 990 GHz spectral range. The dashed red lines mark the position of the detected lines. The 
error bars indicate the median la uncertainty of each spectrum. 



Table 4 

SPIRE/FTS HF and OIT 1 



Fluxes 



Galaxy 


Transition 


forest 


Flux 






(GHz) 


(10" 15 crgcm- 2 s- 1 ) 


UGC 05101 


HF J = 1 - 


1232.476 


-10.4 ± 2.3 


NGC 7130 


HF J = 1 - 


1232.476 


8.6 ± 1.1 




OH+ 1 - 0i 


909.159 


7.6 ± 1.2 




OH+ 1 2 - 0i 


971.805 


9.6 ± 2.4 




OH+ li - 0i 


1033.118 


11.2 ± 1.5 



Note. — Fluxes and la statistical uncertainties for the 
HF and OH + detections. 



Survey (PI: C. D. Wilson) and the Herschel Reference 
Survey (PI: S. Eales) programs and the data were avail- 
able in the Herschel archive. The SPIRE photometer 
consists of three bolometer arrays that observe in three 
spectral bands centered at 250, 350, and 500 fj,m. De- 
pending on the angular size of the galaxies they were 
observed in the small or in the large scan-map mode in 
the three bands. At least two scans in perpendicular 
directions were performed for each galaxy. 

First we processed the raw data using HIPE version 9 
to create the flux calibrated timelines for each bolome- 



ter. The standard HIPE pipeline corrects for instrumen- 
tal effects and attaches the pointing information to the 
timelines. Then the timelines were combined to create 
the maps using Scanamorphos version 18 (Roussel 2012). 
Scanamorphos detects and masks glitches in the data, 
performs the baseline subtraction, and projects the time- 
lines in a spatial grid. 

To measure the nuclear fluxes of the galaxies, first 
we converted the maps produced by Scanamorphos in 
Jybeam -1 to Jyarcsec -1 . Wc adopted the beam areas 
given by the SPIRE Observer's Manual (423, 751, and 
1587arcsec 2 for the 250, 350, and 500 /im bands, respec- 
tively). Then we performed aperture photometry using 
an aperture of diameter 18" in the 250 /xm image and 
30" in the 350 and 500 /im maps and we multiplied the 
measured fluxes by an aperture correction. The values 
of these corrections are 2.33, 1.75, and 3.02 for the 250, 
350, and 500 /im fluxes, respectively. These aperture cor- 
rection factors were calculated from a map of Neptune 
processed following the same reduction steps described 
above. We did not apply any color correction to the 
fluxes because this correction is expected to be small, 
<5%, for black-body emission (SPIRE Observer's Man- 
ual). The measured fluxes are presented in Table 5. 

For most of these galaxies, the nuclear region is only 
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Figure 2. (Continued) 



marginally resolved at the angular resolution of SPIRE 
(18", 25", and 37" for the 250, 350, and 500 /im bands, 
respectively). However, for NGC 4051 and NGC 4151 
the nuclear emission is more extended, so the reported 
nuclear flux for these two galaxies might be overesti- 
mated due to the aperture correction applied. For the 
two more distant galaxies in our sample, UGC 05101 
and NGC 7130, the fluxes in Table 5 correspond to the 
galaxy integrated sub-mm emission, since they are point 
sources in the SPIRE images. For the rest, the fluxes 
should be considered a lower limit to their total sub-mm 



emission. 



3.3. Ground-based CO Observations 



We performed pointed observations of low-J CO 
lines in three galaxies, UGC 05101, NGC 3227, and 
NGC 3982. The observations were conducted from 
2011 October 1st to October 3rd with the IRAM 
30m telescope on Pico Veleta (Spain). We used the 
heterodyne eight mixer receiver (EMIR) using differ- 
ent frequency setups, depending on the target. For 
NGC 3227 and NGC 3982 we used two frequency se- 
tups, one centered at 98.875 GHz for the 3 mm band 
(E0) to encompass the emission of CO J = 1 — and 
13 CO J = 1 — 0, and another one that combines the 
3 mm band (E0) and 1mm band (E2), centered at 
89.375 GHz and 230.538 GHz, respectively, to cover the 



Table 5 

SPIRE Photometry 



Galaxy 


F„(250Ami) 


F„(350 fim) 


F„(500 /W) 




(Jy) 


(Jy) 


(Jy) 



NGC 1056 


2.15 


± 


0.03 


1.09 


± 


0.02 


0.439 


± 


0.012 


UGC 05101 


5.83 


± 


0.11 


2.35 


± 


0.04 


0.785 


± 


0.019 


NGC 3227 


3.14 


± 


0.02 


1.55 


± 


0.01 


0.611 


± 


0.005 


NGC 3982 


1.89 


± 


0.01 


1.24 


± 


0.01 


0.566 


± 


0.009 


NGC 4051 


1.49 


± 


0.02 


0.780 


± 


0.011 


0.343 


± 


0.011 


NGC 4151 


0.592 


± 


0.007 


0.305 


± 


0.005 


0.133 


± 


0.008 


NGC 4388 


2.39 


± 


0.07 


1.33 


± 


0.03 


0.614 


± 


0.019 


IC 3639 


1.56 


± 


0.02 


0.845 


± 


0.011 


0.392 


± 


0.012 


NGC 7130 


6.49 


± 


0.10 


2.92 


± 


0.04 


1.04 


± 


0.02 


NGC 7172 


4.14 


± 


0.06 


1.84 


± 


0.03 


0.711 


± 


0.015 


NGC 7582 


16.2 


± 


0.3 


6.47 


± 


0.09 


2.26 


± 


0.04 



Note. — Nuclear fluxes measured in the SPIRE 250, 
350, and 500 /im images. The aperture size is 18" for 
the 250 fim flux and 30" for the 350 and 500 /J,m fluxes. 
Aperture corrections have been applied (see Section 3.2). 



emission of the CO J = 2 - 1 line. For UGC 05101 the 
receiver was tuned at 98.875 GHz and 230.538 GHz to 
cover the CO J = 1 - and CO J = 2 - 1 lines. The 
13 CO J = 1 — was outside the observed bands for 
UGC 05101. The 16 GHz bandwidth of the EMIR re- 
ceivers were connected to the 195 kHz resolution Fourier 
transform spectrometers (FTS), providing a channel- 
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width spacing of ^0.6 km s" 1 for the 3 mm band and 
^0.3 km s" 1 for the 1mm band, and to the wideband 
line multiple autocorrelator (WILMA) that provides a 
resolution of 2 MHz (5.2 — 6.7 km s -1 for the 3 mm band 
and ^2.6kms _1 at 230.538GHz). The half power beam 
width (HPBW) is 21'.'3 at 115 GHz and 10'.' 7 at 230 GHz. 
Observations were performed in wobbler switching mode, 
with a throw of 120". Weather conditions were good, 
with zenith optical depths between 0.1 and 0.3. System 
temperatures ranged from 100 K to 200 K in the 3 mm 
band and from 250 K to 800 K in the 1mm band. Point- 
ing was checked every two hours and the rms pointing 
error was 4" at 3 mm. Focus was checked at the begin- 
ning of each run and during sunrise. 

The data were reduced using the continuum and 
line analysis single-dish Software (CLASS) package of 
the Grenoble image and line data analysis software 
(GILDAS) 15 . The observed spectra were converted from 
antenna temperature units (T*) to main beam tempera- 
tures (T mb ) using the relation T mh = T* x (F eS /B eS ), 
where F c g is the forward efficiency of the telescope 
(i*eff = 0.95 for the 3 mm band and F c s = 0.92 for the 
1 mm band) and B e g is the main beam efficiency, rang- 
ing from 0.77 to 0.81 at 3 mm, and B e g = 0.59 at 1 mm. 
A second order baseline was removed for each individ- 
ual scan and co-added after inspection. The final spec- 
tra were smoothed to a common velocity resolution of 
20kms~ 1 . In this work we present the results of the CO 
and 13 CO lines. The other molecular species detected 
within the EMIR bands will be presented and analyzed 
in a forthcoming paper. 

In addition, we analyzed the CO J = 3 — 2 transition 
for these galaxies observed with the 15 m James Clerk 
Maxwell Telescope (JCMT) on Mauna Kea (USA) and 
the 10 m Hcinrich Hertz Telescope (HHT) on Mt. Gra- 
ham (USA). These observations are described in detail by 
Wilson et al. (2008), Iono et al. (2009), and Mao et al. 
(2010). 

3.3.1. CO Line Profile Fitting 

The observed low-J CO transitions have broad and 
complex line profiles suggesting that the CO emission is 
produced in several regions unresolved within the tele- 
scope beam. To estimate the CO intensity arising from 
each region we fitted multiple Gaussians to the observed 
line profiles. Three Gaussians reproduce well the ob- 
served profiles in these galaxies. 

First, we derived the FWHM and relative velocity of 
each kinematical component by fitting simultaneously 
the high signal-to-noise (S/N) ratio CO J =1 — and 
CO ,7 = 2 — 1 line profiles. We assumed that the FWHM 
and relative velocity of each kinematical component are 
equal for both transitions, and only the peak intensities 
were allowed to vary. Then, to model the lower S/N ra- 
tio CO J = 3 - 2 and 13 CO J = 1 - lines we used the 
FWHM and velocity values derived before and fitted the 
peak intensities of each component. The best-fit models 
are shown in Figures 3, 4, and 5. In Table 6 we show the 
results of the Gaussian fits. 

3.3.2. Literature CO Data 
15 http://www.iram.fr/IRAMFR/GILDAS 



UGC05101 



0.03 
g 0.02 
.1 0.01 




11000 11500 12000 12500 
v LS r [ km s" 1 ] 

Figure 3. Observed line profiles and best-fit models of the 
CO J = 1 - and CO J = 2 - 1 transitions observed with IRAM 
30 m and the CO J = 3 — 2 transition observed with JCMT 15 m 
for UGC 05101. The spectra are smoothed to a velocity resolution 
of 20kms~ 1 . The narrower spectral coverage of the CO J = 3 — 2 
transition is due to the narrower instantaneous bandwidth of the 
JCMT receivers. 

We obtained from the literature the observed in- 
tensities of the CO ,7=1-0, CO ,7 = 2-1, and 
CO J = 3 — 2 transitions for the remaining galaxies (see 
Table 7). For most of them the CO J = 1 - and 
CO J = 2 — 1 intensities are available. However, the 
CO J = 3 — 2 transition was observed in just three of 
them (see also Section 3.3.1). 

For the three galaxies analyzed in Section 3.3.1, in Ta- 
ble 7 we list the CO fluxes corresponding to the compo- 
nent with the central recessional velocity. We associate 
this component with the compact nuclear CO emission, 
while the other two components are likely due to the 
rotation of the galaxy disk. The disk emission also con- 
tributes to the central component, so the adopted nuclear 
flux might be slightly overestimated. 

4. CO EMISSION MODELING 

In this section we describe the analysis of the CO 
spectral line energy distributions (SLEDs) of six of our 
galaxies with more than three CO lines detected in their 
SPIRE/FTS spectra. First we describe the radiative 
transfer models used to interpret the CO SLEDs. The 
models for individual sources are discussed later. Finally 
we derive the properties of the cold and warm molecular 
gas traced by the CO lines. 

4.1. Radiative Transfer Models 

We used the non local thermodynamic equilib- 
rium (non-LTE) radiative transfer code RADEX 



Sub-millimeter Spectra of Local Active Galaxies 



9 



Table 6 

CO Line Profile Fitting Results 



— 


~ 

"LSR 


r VV xllvl 






pn j — q o 
V_, \J J — o — z 






(kms" 1 ) 


(kms" 1 ) 


(Kkms- 1 ) 


(Kkms- 1 ) 


(Kkms- 1 ) 


(Kkms- 1 ) 


UGC 05101 


11784 


318 ± 32 


9.64 ± 0.38 


19.46 ± 0.64 


7.9 ± 1.6 






11581 


144 ± 8 


2.58 ± 0.17 


7.07 ± 0.29 


3.53 ± 0.72 






12025 


175 ± 5 


4.45 ± 0.21 


11.04 ± 0.35 


6.22 ± 0.87 




NGC 3227 


1110 


206 ± 8 


31.83 ± 0.45 


53.17 ± 0.53 


15.4 ± 2.0 


1.70 ± 0.15 




972 


75 ± 7 


2.97 ± 0.16 


4.92 ± 0.19 


< 2.2 


< 0.2 




1295 


122 ± 4 


7.21 ± 0.27 


17.60 ± 0.31 


< 3.6 


0.53 ± 0.09 


NGC 3982 


1114 


102 ± 9 


7.31 ± 0.26 


8.53 ± 0.33 


4.3 ± 1.4 


0.70 ± 0.10 




1061 


46 ± 5 


3.25 ± 0.12 


4.20 ± 0.15 


< 1.8 


0.28 ± 0.05 




1189 


45 ± 8 


1.87 ± 0.12 


< 0.5 


4.94 ± 0.61 


0.14 ± 0.05 



Note. — The uncertainties of the measured radio local standard of rest (LSR) velocities 
are ~30kms _1 . 



Table 7 

CO Line Fluxes from Ground-based Observations 



Galaxy 


COJ= 1-0 
(10- 16 erg cm" 2 s" 1 ) 


Ref. a 


CO J = 2- 1 
(10- 16 erg cm" 2 s" 1 ) 


Ref. a 


CO J = 3-2 
(TO -16 erg cm" 2 s- 1 ) 


Ref. a 


NGC 1056 


2.5 ± 0.5 


1,1 


10.4 ± 2.1 


1,1 






UGC 05101 


1.8 ± 0.4 


2,1 


9.7 ± 2.0 


2,1 


22.5 ± 6.4 


3,J 


NGC 3227 


6.1 ± 1.2 


2,1 


26.5 ± 5.3 


2.1 


125 ± 30 


4,2,H 


NGC 3982 


1.4 ± 0.3 


2,1 


4.3 ± 0.9 


2.1 


34.9 ± 13.0 


4,2,H 


NGC 4051 


5.6 ± 1.2 


5,N 










NGC 4151 






3.0 ± 1.2 


6.1 






NGC 4388 


1.7 ± 0.4 


7.1 


5.5 ± 1.7 


7.1 






IC 3639 


2.2 ± 0.4 


l.S 


20.9 ± 4.2 


l.S 






NGC 7130 


12.4 ± 2.5 


l.S 


44.8 ± 9.0 


l.S 






NGC 7172 


13.3 ± 2.7 


l.S 


26.1 ± 5.2 


l.S 






NGC 7582 


18.0 ± 3.7 


8,S 


139 ± 30 


8,S 







References. - (1) Albrecht et al. (2007). (2) This work (Section 3.3). (3) Wilson et al. 
(2008). (4) Mao et al. (2010). (5) Vila-Vilaro et al. (1998). (6) Rigopoulou et al. (1997). (7) 
Pappalardo et al. (2012). (8) Aalto et al. (1995). 

Note. - Fluxes and la uncertainties of ground-based observations of the low-J CO lines. 
The fluxes were converted from their original units (Kkms" 1 ) using the following relation: 
F(ergcm" 2 s _1 ) = /(Kkms" 1 ) 9.197 x 10~ 13 iy(GRz)/(irri MB (R(cm)) 2 ), where v is the transition 
frequency, t/mb the main beam efficiency, and R the radius of the telescope. 

(a) Reference for the CO line flux and telescope used. (I) IRAM 30 m; ( J) JCMT 15 m; (H) HHT 
10 m; (N) NRO 45 m; (S) SEST 15 m. 



(van der Tak et al. 2007) to infer the physical conditions 
in the regions where the mid- J CO emission is produced. 
RADEX uses the escape probability approximation to 
solve iteratively for the molecular level populations and 
the intensities of the CO lines. 

We produced a grid of models assuming a spherically 
symmetric homogeneous medium. These models span a 
wide range in kinetic temperature (Tkin = 50 — 3000 K), 
molecular hydrogen density (nn 2 = 10 — 10 5 cm- 3 ), 
and CO column density per unit of line width (Nco I 
Av = 10 12 - 10 17 cm- 2 /kms- 1 ). We adopted the colli- 
sional rate coefficients between CO and H2 of Yang et al. 
(2010). For the background radiation we used the cosmic 
microwave background (CMB) at 2.73 K, the inclusion of 
the local far-IR background does not modify significantly 
the models (see Spinoglio et al. 2012). 

The CO lines are spectroscopically unresolved in the 
FTS spectra of these galaxies, so we used ground-based 
observations of the low-J CO lines to determine the value 
of Av (see Section 4.2 for details). 

The free parameters of the models are the T^n, nH 2 , 



Aco, an d the filling factor The filling factor rep- 
resents the fraction of the beam area occupied by the 
source. When the emission is optically thin (t <C 1), 
in the escape probability approximation, the line ratios 
only depend on Tki n and nn 2 , whereas the absolute fluxes 
depend on the Nco® product. Consequently, it is not 
possible to determine simultaneously Nqo aid the fill- 
ing factor, and since the size of the CO emitting region 
(i.e., the filling factor) in our galaxies is not accurately 
known, we can only estimate the beam-averaged column 
density, <A r C o>= : N C o$- 

We estimated the values and uncertainties on the 
model parameters using a Monte Carlo method. From 
the fluxes and their associated uncertainties listed in Ta- 
bles 3 and 7 we constructed 500 CO SLED simulations 
for each galaxy that were fitted using a x 2 minimization 
algorithm. The parameter values and la uncertainties 
reported in Table 8 are those obtained from the param- 
eter distributions in the Monte Carlo simulations. The 
best-fit models of the six analyzed galaxies are shown in 
Figure 6. 
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Figure 4. Same as Figure 3, but for NGC 3227. For this galaxy 
the CO J = 3 — 2 transition was observed with HHT 10 m. The 
narrower spectral coverage of the CO J = 3 — 2 transition is due to 
the narrower instantaneous bandwidth of the HHT receivers. 
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Figure 5. Same as Figure 3, but for NGC 3982. For this galaxy 
the CO J = 3 — 2 transition was observed with HHT 10 m. 



From the derived column densities we calculated the 
molecular gas masses using the following expression: 



mol 



\i <Nco> xco 1 £ld 2 mH 2 



(i) 



where /i = 1.4 is the mean molecular weight per H 
atom, < iVco > is the beam-averaged CO column density, 
x co = 3 x 10~ 4 is the assumed CO abundance with re- 
spect to H2 in dense molecular clouds (Lacy et al. 1994; 
Hernandez & Tan 2011), Q is the adopted SPIRE SLW 
beam size (1.66 x 10 -8 sr), d is the distance to the galaxy, 
and mn 2 is the mass of an H 2 molecule. In Table 9 we 
present the calculated molecular gas masses. 

4.2. Analysis of Individual Sources 
4.2.1. UGC 05101 



^• 95 L ) 



UGC 05101 is an IR bright galaxy (L m = 10 1 
located at 163 Mpc (Sanders et al. 2003). It has a pe- 
culiar morphology that suggests a recent interaction 
(Sanders et al. 1988). Its nuclear activity is classified 
as LINER from optical spectroscopy (Yuan et al. 2010). 
However the detection of high- ionization emission lines in 
the mid-IR ([Ne V] and [O IV]; Armus et al. 2007), the 
intense FeKa emission line at 6.4 keV (Imanishi et al. 
2003), and its high hard X-ray luminosity (£ 2 _iokeV = 
7.5 x 10 42 ergs _1 ; Imanishi et al. 2003) indicate the pres- 



ence of a buried AGN with a luminosity comparable to 
that of a Seyfert galaxy. 

In the SPIRE/FTS spectrum of this galaxy we detect 
five unblended CO lines (see Table 3). An emission line 
is detected at the frequency of the COJ=10 — 9 tran- 
sition, but it is likely blended with the 0-H2O 3i2~22i 
1153 GHz line since water and CO emission lines have 
comparable fluxes in this galaxy. With the spectral res- 
olution of SPIRE/FTS it is not possible to deblend the 
two emission lines; therefore we excluded it from the ra- 
diative transfer analysis. 

In addition to the SPIRE/FTS data, we used ground- 
based observations of the low- J CO transitions to con- 
strain the radiative models. Interferometric observa- 
tions of the CO J = 1 — line show that the CO emis- 
sion is produced in a compact region (<lkpc) plus a 
fainter extended (2.2 kpc) rotating disk (Genzel et al. 
1998; Wilson et al. 2008). 

The observed FWHM of the CO J = 3 - 2 line is 
650 km s" 1 (Iono et al. 2009). However, the SPIRE/FTS 
mid- J CO lines are spectroscopically unresolved, that is, 
the width of mid- J lines is <500kms _1 . The line profile 
fitting of the CO J = 1 - 0, J = 2 - 1, and J = 3 - 2 
transitions (Section 3.3.1) shows that the width of the 
middle component is 320 km s" 1 , so we assumed that it 
is also the width of the mid- J CO lines. Consequently, 



Sub-millimeter Spectra of Local Active Galaxies 



11 



Table 8 

RADEX Models 



Galaxy 


Av a 


logT kin b 


log 




log < N CO > d 




e 




(kms- 1 ) 


(K) 


(cm 3 ) 


(cm" 2 ) 


(10 15 ergcm 2 s 


UGC 05101 


320 


2.9 ± 0.3 


3.7 


± 0.4 


15.2 ± 0.1 


102 ± 


li 


NGC 3227 


200 


2.8 ± 0.1 


3.3 


± 0.1 


15.8 ± 0.1 


190 ± 


10 


NGC 3982 


100 


2.0 ± 0.1 


1.5 


± 0.4 


15.2 ± 0.2 


39 ± 


4 


NGC 4388 


125 


2.5 ± 0.2 


3.9 


± 0.2 


15.4 ± 0.1 


104 ± 


12 


NGC 7130 


90 


2.8 ± 0.1 


3.2 


± 0.1 


16.1 ± 0.1 


240 ± 


10 


NGC 7582 


200 


2.6 ± 0.2 


3.6 


± 0.7 


16.4 ± 0.2 


540 ± 


30 



Note. — Best-fit parameters of the RADEX models. Line width assumed 
for the model. Kinetic temperature. Molecular hydrogen density. Beam- 
averaged CO column density. The beam FWHM is 30". ' e) Integrated flux of the 
CO lines from J up = 4 to J up = 12 of the best-fit model. 



we used the CO J = 1 - 0, J = 2-1, and J = 3-2 fluxes 
of this middle component in the CO SLED fit presented 
in Figure 6. 

4.2.2. NGC 3227 

NGC 3227 is a nearby (14.4 Mpc) Seyfert 1.5 galaxy 
(Ho et al. 1997). It is interacting with the dwarf ellip- 
tical galaxy NGC 3226. The hard X-ray luminosity of 
the AGN is 1.9xl0 41 ergs" 1 (Gondoin et al. 2003).' In 
addition, recent (^40 Myr) intense star-formation is ob- 
served in the nucleus. The luminosity produced by these 
young stars represents about half of the bolometric lumi- 
nosity of the galaxy (Davies et al. 2006). The Spitzer/ 
IRS mid-IR spectrum of this galaxy shows the [Ne V] 
and [O IV] high-ionization emission lines confirming the 
presence of an AGN (Deo et al. 2007; Dasyra et al. 2008; 
Weaver et al. 2010). 

Intcrferometric observations of this galaxy 
(Schinnerer et al. 2000) show that part of the 
CO J = 1 — and J = 2 — 1 emission is concentrated in 
a compact region of ~8" of diameter (~500pc). 

In the FTS spectrum of NGC 3227 we detected the 
CO lines from J up = 4 to J up = 11. These CO emis- 
sion lines are not resolved at the SPIRE/FTS spectral 
resolution; thus we assumed the line width of the middle 
component of the CO J = 1 — and J = 2 — 1 transitions, 
200kms- 1 (sce Section 3.3.1). 

4.2.3. NGC 3982 

NGC 3982 is a Seyfert 1.9 (Ho et al. 1997) located 
at 20.6 Mpc. The hard X-ray luminosity of this galaxy 
is <2.5xl0 39 ergs _1 and the FeKa feature at 6.4keV 
is not detected (Guainazzi et al. 2005). This suggests 
that NGC 3982 hosts a relatively low luminosity AGN. 
Moreover, the luminosities of the high-ionization mid-IR 
[O IV] and [Ne V] emission lines are about 2xl0 39 ergs _1 
(Tommasin et al. 2010), low compared to those of other 
Seyfert galaxies (see Figure 2 of Pereira-Santaella et al. 
2010) and confirming that the AGN in NGC 3982 is not 
very luminous. 

Four CO lines, from J up = 4 to J up = 7, are detected in 
the FTS spectrum of this galaxy, all of them in the SLW 
range. In the modeling we included ground-based obser- 
vations of the CO J =1 — 0, J = 2 — 1, and J = 3 — 2 
lines (see Section 3.3). The ground- based data shows 
three kinematical components (see Figure 5 and Section 
3.3.1). The width of the middle component is 100 kms . 



The other two components are probably produced by a 
rotating disk. The difference in the relative intensities 
of the rcdshifted and blueshifted disk components from 
the J = 2 — 1 to the J = 3 — 2 transition are presum- 
ably due to pointing uncertainty and the smaller beam 
size of the CO J = 2 — 1 observation. Mapping observa- 
tions of these transitions would be needed to confirm this 
interpretation. 

The best-fit model for NGC 3982 has a higher density 
and lower temperature than the rest of the galaxies (see 
Table 8 and Figure 6). No CO lines with J up > 7 are 
detected, so our fit could be biased toward low kinetic 
temperatures. However, the CO J = 4 — 3/CO J = 
7 — 6 ratio in NGC 3982 is more than a factor of two 
higher than in the other galaxies in the sample. This 
implies a lower kinetic temperature for a fixed density 
and, therefore, the estimated lower temperature is likely 
real and not a consequence of the limited sensitivity of 
the SPIRE/FTS spectrum. 

4.2.4. NGC 4388 

NGC 4388 is a nearby edge-on spiral galaxy located 
at 18.2 Mpc. A weak broad Ha component has been re- 
ported by Shields & Filippenko (1996) and it is classified 
as Seyfert 1.9 by Ho et al. (1997). However Ruiz et al. 
(1994) did not find evidence for broad components in the 
near-IR Pa/3 line. A hidden broad line region (HBLR) is 
detected in optical polarized light (Young et al. 1996). 
Both the mid-IR [Ne V] and [O IV] high-ionization 
emission lines are detected in the Spitzer /IRS spec- 
trum of this galaxy (Tommasin et al. 2010). The ab- 
sorption corrected hard X-ray luminosity of the AGN is 
3.0xl0 41 ergs- 1 (Cappi et al. 2006). A radio jet is ob- 
served in this galaxy (Falcke et al. 1998). 

We detect six CO lines from J up = 4 to 9 in the SLW 
SPIRE/FTS spectrum. We obtained ground-based mea- 
surements of the CO J = 1 — and J = 2 — 1 fluxes 
from the literature (Pappalardo et al. 2012, see Sec- 
tion 3.3). For the line width we assumed 125kms _1 
(Vila-Vilaro et al. 1998). The best-fit model is shown 
in Figure 6. 

4.2.5. NGC 7130 

NGC 7130 is a luminous IR galaxy (L m = 10 1135 L ) 
located at 66 Mpc (Sanders et al. 2003). Its nuclear ac- 
tivity is classified as Seyfert 2 (Yuan et al. 2010). It 
is a face-on spiral galaxy with a disturbed morphol- 
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Figure 6. CO SLED (blue circles) and best-fit model (red line) for our galaxies. The shaded region is the la uncertainty in the best-fit 
model calculated from the Monte Carlo simulations (see Section 4.1). 
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ogy. Faint tidal tails arc observed, although there are 
no nearby companions. It hosts a Compton-thick AGN 
with an estimated intrinsic hard X-ray luminosity of 
8.3 xlO 40 ergs" 1 (Lcvenson ct al. 2005). Spitzer/IRS ob- 
servations show the detection of the [Ne V] and [O IV] 
high-ionization emission lines (Tommasin et al. 2010). In 
addition to the powerful AGN, a strong burst of star- 
formation is taking place in the nuclear region of this 
galaxy (Gonzalez Delgado et al. 2001; Diaz-Santos et al. 
2010). 

In the SPIRE/FTS spectrum of NGC 7130 we mea- 
sured nine CO lines from J up = 4 to 12. We fitted 
the lines observed with SPIRE/FTS together with the 
CO J = 1 — and CO J = 2 — 1 fluxes obtained from the 
literature (Albrecht et al. 2007). The CO line width is 
-^Okins- 1 (Albrecht et al. 2007). 

The CO J = 10 — 9 line was excluded from the fit be- 
cause it could be blended with the H2O312-221 line. 
However the best-fit model explains the observed flux of 
the CO J = 10 — 9 line, and therefore the contribution of 
the water line is likely small (see Figure 6). 

4.2.6. NGC 7582 

NGC 7582 is a Seyfert galaxy at a distance of 
20.6 Mpc. It is a member of the interacting Grus 
Quartet (de Vaucouleurs 1975). The nuclear activity 
is classified as Type 2 because no broad lines are de- 
tected in the optical (Heisler et al. 1997). However, 
Reunanen et al. (2003) reported a broad component in 
the ncar-IR Br7 line. The hard X-ray luminosity of the 
AGN is 2.0xl0 41 ergs" 1 (Piconcelli et al. 2007). The nu- 
clear region of NGC 7582 harbors young star-forming 
regions that dominate the ionization of the nuclear gas 
(Riffcl et al. 2009). In the mid-IR Spitzer/IRS spectrum 
of NGC 7582 the [Ne V] and [O IV] high-ionization emis- 
sion lines are detected (Tommasin et al. 2010). 

The CO lines from J up = 4 to 12 are detected in the 
SPIRE/FTS data. We added to the CO SLED ground- 
based observations of the CO J = 1 — and J = 2 — 
1 lines (Aalto et al. 1995). For the line width we used 
-200kms _1 (Baan et al. 2008). 

The nuclear far-IR emission of this galaxy is only par- 
tially resolved at the resolution of Herschel, so the effects 
of the different SPIRE/FTS beam sizes are minimized. 
Therefore, we fitted all the CO lines simultaneously (Fig- 
ure 6). 

4.3. Cold Molecular Gas 

From Figure 6 we see that the best-fit models for the 
warm molecular gas components (1km ^500 K; see Sec- 
tion 4.4) underpredict the flux of the CO J = 1 — line. 
The predicted fluxes are 10-25 % of the observed values. 
This CO J = 1 — excess is explained if a cold gas com- 
ponent is also present in these galaxies. 

The beam-averaged H2 column density of the cold 
component can be estimated using the CO-to-H2 conver- 
sion factor V H2 (cixr 2 ) = 0.5 x 10 20 / C oj=i-o (Kkms" 1 ) 
revised for active galaxies (Downcs & Solomon 1998). 
Substituting in Equation 1 we obtain: 
M mol (M ) ~ 7 x lO^Fcojrf-oCergcm-^-^x 
(d(Mpc)) 2 . We subtracted the contribution of the warm 
component to the CO J = 1 — line for the six galaxies 
modeled. For the rest of the galaxies this contribution 



to the CO J = 1 — emission introduces an extra ~25 % 
uncertainty in the cold molecular mass. The calculated 
masses are listed in Table 9; they range from 10 7 6 to 
1O 95 M . 

4.4. Warm Molecular Gas 

From the RADEX modeling of the low- and mid- J CO 
lines, we find that the CO SLEDs of five out of the six 
galaxies analyzed (all except NGC 3982) can be fitted 
with a warm gas component with comparable physical 
properties, Hh 2 ~10 3 2 -10 3 - 9 cm" 3 and T kin -300-800 K 
(see Table 8). 

Similar analysis of the mid-J CO SLED have 
been presented for Arp 220, M 82, and NGC 1068 
(Rangwala et al. 2011; Kamenetzky et al. 2012; 
Spinoglio et al. 2012). In these galaxies two or 
more CO components are fitted to the SLED; therefore 
we compare our results with their warm gas component. 
The gas density and temperature in our galaxies are 
close to those of the starburst M 82 (tih 2 = 10 3 ' 4 cm~ 3 
and T^n = 500 K). In the ultraluminous IR galaxy 
(ULIRG) Arp 220 the warm molecular gas has higher 
temperature, T^m = 1300 K, and slightly lower density, 
«h 2 = 10 3 ' 2 cm~ 3 , although within the uncertainties, 
both parameters are comparable to the values of our 
galaxies. In contrast, the warm molecular gas traced 
by the mid-J CO lines in NGC 1068 has a lower 
temperature, 90 K, and higher density, 10 46 cm~ 3 . 

We use Equation 1 to calculate the molecular mass. 
The resulting masses are listed in Table 9. They range 
from 10 6 to 10 8 M Q . These masses are lower than 
the warm molecular component observed in Arp 220, 
10 8 ' 7 M Q , but resemble those observed in M 82 and 
NGC 1068, 1O 62 M and 10 74 M o , respectively. For 
these six galaxies the mass of the cold molecular gas (Sec- 
tion 4.3) is between 40 and 120 times higher than that of 
the warm component, whereas in Arp 220 and M 82 the 
cold to warm molecular gas mass ratio is smaller, — 10. 

Warm molecular gas is also traced by the mid-IR 
rotational H2 lines. The gas temperatures derived 
from the lowest H2 rotational transitions are 100-1000 K 
(Rigopoulou et al. 1999; Roussel et al. 2007) and are 
similar to the temperature range obtained from our ra- 
diative transfer analysis. Therefore, the mid-J CO lines 
and the rotational H2 lines might trace the same warm 
molecular clouds. 

The rotational H2 — S(l) line at 17.03 /xm is ex- 
pected to account for 15-80% of the total H 2 emis- 
sion arising from the warm molecular gas (based on the 
Shaw et al. 2005 models) depending on the gas tem- 
perature and density. Indeed, the H2 S(l) is usually 
the brightest H2 line in the mid-IR Spitzer/IRS spec- 
tra of galaxies (e.g., Roussel et al. 2007). For most of 
the galaxies in our sample H2 S(l) fluxes were presented 
by Tommasin et al. (2008, 2010) 16 . To correct for dust 
extinction first we estimated using the relation A-&/ 
T9.7 = 1.48 (McClure 2009), where r 9 . 7 is the optical 

16 For UGC 05101 the H 2 S(l) flux was published by 
Farrah et al. (2007). For NGC 3227 and NGC 4151, S. Tom- 
masin provided us with the H2 S(l) fluxes measured from the high 
spectral resolution Spitzer/IRS spectra following the method de- 
scribed in Tommasin et al. (2010). The H 2 S(l) fluxes are (26.8 ± 
0.2)xl0" 14 and (5.44 ± 0.82) X 10" 14 erg cm- 2 s" 1 for NGC 3227 
and NGC 4151, respectively. 



14 



Pereira-Santaella et al. 



Table 9 

Molecular Gas Properties 



Galaxy 


log M w am 


Inn- T & 


Inn- T c 
IT \ 


^CO/ iu warm 


Inn M e 

l°g iw cold 


log iWQ J 


<r>~ Q 
XQ" 




\ M o) 


(J \ 


( T / A/f ^ 

(Lq/Mq) 


( M <3) 


t A/f *\ 




NGC 1056 






6.1 




8.3 


4.6 


5.6xl0~ 5 


UGC 05101 


7.7 


7.9 


8.0 


1.4 


9.5 


6.4 


1.5xl0" 4 


NGC 3227 


6.2 


6.1 


6.3 


0.6 


8.0 


4.7 


4.1X10" 4 


NGC 3982 


5.9 


5.7 


5.6 


0.6 


7.6 


1.1 


1.3xl0" 4 


NGC 4051 






5.8 




8.0 


1.0 


2.5xl0~ 5 


NGC 4151 






5.4 






3.8 




NGC 4388 


6.0 


6.0 


6.3 


1.0 


7.6 


1.7 


8.4xl0~ 4 


IC 3639 






6.4 




8.5 


1.8 


4.8xl0~ 5 


NGC 7130 


7.8 


7.5 


7.1 


0.4 


9.6 


6.0 


1.7xl0~ 4 


NGC 7172 






6.8 




9.0 


5.8 


3.6xl0~ 4 


NGC 7582 


7.0 


6.9 


6.7 


0.6 


8.7 


5.1 


1.8xl0" 4 



Note. — For the warm and cold molecular masses (M warm and M co id) listed in this table we 
assumed Af mo i = 1.4Mh 2 to account for He mass. 

^ Mass of the warm molecular gas traced by the mid- J CO lines (see Section 4.4). W Integrated 
luminosity of the CO lines from J up = 4 to J up = 12 of the best-fit RADEX model (see also Table 
8). ' c ' Extinction corrected luminosity of the H2 - S(l) line at 17.03 /im. Warm molecular 
gas cooling rate per unit gas mass due to the CO transitions from J up = 4 to J up = 12. ^ Mass of 
the cold molecular gas derived from the CO J = 1 — line flux (see Section 4.3). ( - 9 ' Neutral carbon 
mass in the cold molecular clouds (Section 5). Neutral carbon abundance with respect to H2. 



depth of the 9.7 /im silicate absorption feature. The 79.7 
values of our galaxies were published by Wu et al. (2009). 
Then, assuming A 17lim ~ 0.4 x A K (McClure 2009), we 
corrected the H 2 S(l) fluxes. In Table 9 we compare 
the extinction corrected luminosity of this line with the 
integrated luminosity of the mid- J CO lines (J up = 4 
to J up = 12, those in the spectral range observed by 
SPIRE/FTS). The total H 2 luminosity is 0.1-0.8dex 
higher than the H2 S(l) luminosity depending on the 
exact gas conditions. So, in general, the H2 cooling is 
higher but comparable to the mid- J CO cooling in these 
galaxies. 

Table 9 also shows the ratio between the integrated 
luminosity of the mid- J CO lines and the molecular gas 
mass (Lco/A^moi)- This ratio is between 0.4 and 1.3 Lq/ 
Mq in our galaxies, similar to the CO cooling ratio mea- 
sured in the warm CO component of Arp 220 (0.4 Lq/ 
Mq Rangwala et al. 2011). For the gas conditions de- 
rived from the modeling, the expected H2 cooling ratio is 
~0.5-10L Q /Af© (Le Bourlot et al. 1999). These values 
agree with the observed rotational H2 cooling; thus this 
supports the idea that the mid-IR rotational H2 lines and 
the mid- J CO lines originate in the same warm molecu- 
lar gas. This agreement also indicates that the CO abun- 
dance (xco = 3 x 10 -4 ) used to calculate the M mo \ is 
reasonable for the warm molecular gas in these active 
galaxies. 



4.5. Warm CO Heating Source 

Our radiative transfer modeling shows that warm 
molecular gas contributes significantly to the CO lumi- 
nosity in these Seyfert galaxies. To determine the heat- 
ing source of the warm gas in Figure 7 wc compare the 
mid- J CO luminosity with the [O IV]25.9^m luminosity, 
a proxy of the AGN luminosity (Melendez et al. 2008), 



and the IR luminosity 17 , which is correlated with the 
star- formation rate (Kennicutt 1998). There is no clear 
correlation between the AGN luminosity and the CO lu- 
minosity (left panel). For .Lroiv] ~ 5 x 10 Lq, the CO 
luminosity differs two orders of magnitude from galaxy 
to galaxy. On the other hand, the Lqo seems consistent 
with a Lco/Liji ratio about 10~ 4 . That is, the mid- J 
CO emission is more likely related to the star-formation 
activity than to the AGN luminosity in our sample of 
Seyfert galaxies. 
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Figure 7. Mid- J CO luminosity vs. [O IV]25.9^tm luminosity 
(left panel) and vs. IR luminosity (right panel). The [O IV]25.9/im 
luminosity is a proxy of the AGN luminosity in Seyfert galaxies 
whereas the IR luminosity traces the star-formation rate. The 
dashed lines in the right panel show the ico/^lR ratio. The 
[O IV]25.9^tm fluxes are from the Pereira-Santaella et al. (2010) 
compilation. 

PDRs are expected in star-forming regions where UV 
photons from young stars heat the ISM; therefore we 

17 The values in Table 1 are calculated from IRAS fluxes. 
IRAS beams are much larger than the SPIRE/FTS beam, thus 
for the galaxies with the larger apparent sizes (see Figure 1) wc 
consider this Ltr as an upper limit. 
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Figure 8. CO SLEDs normalized to the integrated CO luminosity 



(from J up = 1 to J u 



12). We compare the average CO SLED 



of our Seyfert galaxies with those of Mrk 231 (van der Werf et al. 

2010) , NGC 6240 (Meijerink et al. 2013), Arp 220 (Rangwala et al. 

2011) , and M 82 (Kamenetzky et al. 2012). 



compared the observed CO SLEDs with the PDR models 
of Wolfire et al. (2010). However, we find that no single 
PDR model reproduces the observed CO SLED. In gen- 
eral the best-fit PDR models underpredict the fluxes of 
the CO lines with J up > 8 lines. Conversely, XDR mod- 
els (Meijerink et al. 2006) reproduce better the higher J 
CO lines and slightly underpredict those with J up < 6. 
Consequently we cannot rule out an XDR contribution 
to the higher J CO emission in our galaxies. 

In Figure 8 we compare the average CO SLED of our 
Seyfert galaxies with those published for other galaxies. 
The luminosity normalized CO SLEDs arc visibly very 



similar up to J U r 



11. Only for J up > 11, the CO 



emission in Mrk 231 and NGC 6240 is approximately a 
factor of 3 brighter than in the rest of galaxies. For the 
Seyfert 1 Mrk 231 an XDR is proposed to explain the CO 
J up > 10 emission and two PDR to model the lower- J 
CO emission (van der Werf et al. 2010), whereas for the 
merger galaxy NGC 6240 shock excitation can explain its 
complete CO SLED (Meijerink et al. 2013). For Arp 220 
PDR, XDR, and cosmic rays heating mechanisms arc 
ruled out and mechanical energy from supernovae is sug- 
gested as the heating source (Rangwala et al. 2011). For 
M 82 a combination of PDR and shocks is required to 
explain the observed CO emission (Kamenetzky et al. 
2012). That is, the mid-J CO emission of galaxies, at 
least up to J up ~ 11, seems to be produced by warm 
molecular gas with similar average physical conditions, 
although it can be heated by different combinations of 
mechanisms (PDR, XDR, cosmic rays, and shocks). 

5. [C I] 492 AND 809 GHZ EMISSION 

PDR models predict that the two [C I] lines at 492 
and 809 GHz are produced in the transition region from 
C + to CO that occurs in a relatively thin layer within 
the molecular clouds (Kaufman et al. 1999; Bolatto et al. 
1999). However some observational evidence suggests 
that neutral carbon and CO can coexist in the same 
volume (see Papadopoulos et al. 2004 and references 
therein). The latter can be achieved due to non chemi- 
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logT kin [K] 



Figure 9. The dashed black lines are the [Cl]809GHz/ 
[Cl]492GHz ratio (fluxes in ergem" 2 s" 1 sr" 1 ) as a function of 
the H2 density (njj 2 ) and kinetic temperature (T kin ) in the opti- 
cally thin limit. The solid green lines are the ratios measured in 
our sample of galaxies. The green circles represent the molecu- 
lar gas conditions derived from the mid-J CO lines (see Section 
4). UGC 05101 is not included in this figure because only the 
[C I] 809 GHz line is detected in this galaxy. 



cal equilibrium between C° and CO (Papadopoulos et al. 
2004; Glover et al. 2010) or turbulent diffusion of the 
neutral carbon layer into the interiors of the molec- 
ular clouds (Xie et al. 1995). Moreover, cosmic rays 
and X-ray radiation also increase the C°/CO abundance 
ratio throughout the clouds (Papadopoulos et al. 2004; 
Meijerink et al. 2006). 

For five of the galaxies in our sample we detect both 
[C I] fine structure emission lines at 492 and 809 GHz. 
For the rest we detect only the [C I] line at 809 GHz (see 
Table 3). This is in part because the 492 GHz emission 
line lies in a noisy part of the SLW spectra, but also 
because it seems to be weaker than the 809 GHz line (at 
least in the galaxies with both [C I] lines detected). 

To investigate the origin of the [C I] emission in these 
galaxies we used RADEX to calculate the fluxes of the 
[C I] 492 GHz and 809 GHz lines for wide range of physical 
conditions (n, H2 =10-10 8 cm" 3 , T=10-1000K, and N c / 
Aw=10 12 -10 18 cm- 2 /(kms- 1 )). We considered C° col- 
lisions with H2 using the collisional rate coefficients of 
Schroder et al. (1991). 

In Figure 9 wc plot the expected [C I] 809 GHz/ 
[C I]492GHz ratio as a function of the molecular hydro- 
gen density and kinetic temperature for optically thin 
emission. The solid lines represent the observed ratios 
for the five galaxies with both lines detected. They 
span a narrow range in our sample, from 1.2 to 2.2. 
These [C I]809GHz/[C I]492GHz ratios indicate a ki- 
netic temperature about 19-27 K if the gas density is 
higher than the critical densities of the neutral carbon 
lines (tjh 2 > 10 4 cm" 3 ). For lower gas densities, the tem- 
perature is not well constrained since this ratio would be 
compatible with that of higher temperature (Xki n >50K) 
diffuse clouds (lower right corner in Figure 9). 

To test if the [C I] emission arises from the warm molec- 
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ular gas component traced by the mid-J CO lines, we 
plotted the derived physical conditions of the warm gas 
(see Section 4) in Figure 9. Under these conditions the 
expected [C I]809GHz/[C I]492GHz ratio is 6-8, more 
than three times the observed ratio, so the bulk of the 
neutral carbon emission is not likely to be produced in 
the warm molecular gas component. 

Instead, the [C I] emission can originate in the cold gas 
component traced by the CO J = 1 — line (see Section 
4.3). Actually there is a good correlation between the C° 
and CO J = 1 — intensity in our Galaxy (Ikeda et al. 
2002; Ojha et al. 2001) and nearby galaxies (Wilson 
1997). However, the ratio between the CO and [C I] emis- 
sions depends on the specific physical conditions (UV 
flux, cosmic ray rate, turbulence, etc.) in the molecular 
clouds (Israel & Baas 2002; Papadopoulos et al. 2004). 

We used the RADEX grid of models to calculate the C° 
mass and column density in these cold molecular clouds. 
We assumed that the neutral carbon emission is opti- 
cally thin and thermalized (nn, > 10 4 cm -2 ). From 
the [C I]809 GHz/[C I]492 GHz ratio we estimated the gas 
temperature (see Figure 9). When the [C I]492GHz line 
is not detected we assumed T^n = 23 K (the median 
for the galaxies with both [C I] lines detected). Then 
we compared the observed [C I] 809 GHz luminosity with 
the predicted value to estimate the C° column densities 
and masses. The derived C° masses range from 10 3 8 - 
10 6 - 3 M© (Table 9). In Section 4.3 we calculated the mass 
of the cold molecular component; therefore using that the 
relative weight of C° and H2 is 6, we obtain that the C° 
abundance relative to H 2 in these galaxies is ~3xl0 -5 - 
8xl0 -4 . Consequently, the C°/CO abundance ratio 
ranges from 0.1 to 2 (assuming xqo = 3 x 10~ 4 ). These 
C°/CO abundance ratios are comparable to those mea- 
sured in quiescent nearby galaxies (Israel & Baas 2002) 
and starbursts (Rangwala et al. 2011; Kamenetzky et al. 
2012). However, it should be noted that these C° abun- 
dances rely on the CO-to-H2 conversion factor that can 
vary up to a factor of ~4 from galaxy to galaxy (e.g., 
Downes & Solomon 1998). 

The C° abundance is expected to be enhanced in XDRs 
with respect to that in PDRs (Meijerink et al. 2006). 
For a constant gas density Meijerink et al. (2006) mod- 
els predict higher [C I]492 GHz/CO J = 1 - intensity 
ratios in XDRs than in PDRs. In these PDR models the 
[C I]492 GHz/CO J = 1 - ratio ranges between 20 and 
170, decreasing with increasing gas densities, whereas 
in XDR models this ratio is between 50 and 10 5 . In 
most of our galaxies the [C I]492 GHz/CO J = 1 - ra- 
tios are between 15 and 40, so they are compatible with 
PDR excitation. Only in NGC 4388 the [C I]492GHz/ 
CO J = 1 — ratio is higher, 85±20. This ratio can be 
explained by PDR models, but also by a high-density 
(n = 10 5 cm~ 3 ) XDR model. Therefore, in most of our 
sample the C° emission is consistent with that expected 
from PDRs, although for one galaxy XDR excitation is 
possible. 

6. HYDROGEN FLUORIDE 

The HF J = 1 - transition at 1232 GHz has been 
recently detected in the SPIRE/FTS spectra of several 
local galaxies: the ULIRGs Mrk 231 (van der Werf et al. 
2010) and Arp 220 (Rangwala et al. 2011); the star- 



burst M 82 (Kamenetzky et al. 2012); and the Seyfert 
2 NGC 1068 (Spinoglio et al. 2012). This transition has 
also been observed in the Clovcrleaf quasar at z = 2.46 
(Monje et al. 2011). In our sample, we detect this line in 
UGC 05101 in absorption, and in NGC 7130 in emission 
(see Table 4). 

Hydrogen fluoride molecules are rapidly formed by the 
reaction of F with molecular hydrogen when the latter 
becomes abundant, so it is a good proxy for molecular 
gas. 

Most of the fluorine in molecular clouds is present as 
hydrogen fluoride, and in diffuse clouds HF can be several 
times more abundant than CO (Ncufeld et al. 2005). 

Owing to the high critical density of the HF J = 1 — 
transition (n H2 ~10 10 cm -3 at 50K; Neufeld et al. 2005), 
it is expected to be observed in absorption. However, IR 
radiative pumping 18 can populate levels with J > 1, and 
thus it would be possible to observe the J = 1 — transi- 
tion in emission (Neufeld et al. 2005). Alternatively, HF 
can be formed in excited states that lead to the emission 
of the J = 1— line. This mechanism, chemical pumping, 
is possible in strong UV radiation fields where HF pho- 
todissociation and formation rates are high (van der Tak 
2012). In addition, since HF has a large dipolc moment, 
collisions with electrons can excite its rotational levels 
(van der Tak et al. 2012). 

6.1. UGC 05101 

We detect the HF J = 1 — transition in absorption 
in UGC 05101. Assuming that this transition is opti- 
cally thin, we calculated the HF column density using 
the following equation (see Draine 2011): 

»r/ -2n n 00 „ m5 W *( cm ) Kj (GHz)) 4 g; 

iV((cm J = 9.33 x 10 - — -. _ . , (2) 

A u i(s ) g u 

where iV; is the column density of the lower energy level, 
A u i the Einstein coefficient, W\ the equivalent width of 
the line, and g u and gi the degeneracy of the upper and 
lower levels, respectively. For the HF J = 1 — transi- 
tion, A w = 2.42 x 10~ 2 s _1 , gi = 3, and go = 1, and P\o 
= 1232.476 GHz (from the Leiden Atomic and Molecular 
Database [LAMDA]; Schoier et al. 2005). In UGC 05101, 
the continuum at the frequency of the HF J = 1 — line 
is 5.4 ± 0.2 Jy, so W\ = (3.8 ±0.9) x lO^cm. Substitut- 
ing in Equation 2 we obtain N a = (1.1±0.3) x 10 14 cm~ 2 . 
This column density should be considered a lower limit 
because: (1) the measured equivalent width is a lower 
limit, since we assumed that all the observed continuum 
at 1232 GHz illuminates the molecular clouds that pro- 
duce the HF absorption and this might not be true; (2) 
we can only observe the HF in molecular clouds in front 
of the far-IR continuum source in our line of sight; and 
(3) we do not consider HF molecules in excited energy 
levels, although most of the HF molecules are expected 
to be in the ground energy level (see Neufeld et al. 2005). 



18 Ground-state HF can be excited by 2.5 (im near-IR photons to 
its first vibrational level, that then can decay to excited rotational 
levels of the ground vibrational state. 
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6.2. NGC 7130 

In NGC 7130 the HF J = 1 -0 transition is detected in 
emission. In this section we discuss three possible exci- 
tation mechanisms to explain the HF J = 1 — emission 
line in this galaxy. 

First we calculated the column density of the upper 
level assuming that the HF emission is optically thin. 
This value is independent of the excitation mechanism. 
Using RADEX, we verified that the optically thin ap- 
proximation is valid for a wide range of physical con- 
ditions (T kin = 5 - 250 K, tih 2 = 10 - 10 10 cnr 3 , and 
A HF = 10 8 - 10 13 cm- 2 for Av = ^kms" 1 ). Therefore 
we can use following relation to calculate the column 
density of the J = 1 level: 



<N U >= 



4tt R 



Q hv„.iA 



(3) 



where < N u > is the beam averaged column density of 
the upper transition level, fi the beam solid angle, F u i 
the line flux in erg cm -2 s _1 , and h the Planck constant. 
For the HF ,7 = 1 — line the beam FWHM is 18", thus 
<Ni >=(9.1 ± 1.2)xl0 10 cm~ 2 . This column density is a 
lower limit to the total HF column density since we have 
only considered the molecules in the J = 1 rotational 
level. 

6.2.1. IR Pumping 

Carroll & Goldsmith (1981) show that IR pumping can 
efficiently increase the population of ground state rota- 
tional levels when 



> 



.4 



J+i,J 



e hv/kT _ l Av 



(4) 



where / is a factor to account for the geometric dilution 
and dust emissivity, v the frequency of the vibrational 
transition, T is the radiation field temperature, and 
Aj + \ t j and A u>v -i are the Einstein coefficients for the ro- 
tational and vibrational transitions, respectively. For the 
lowest vibrational transition A = 63.2 s _1 and hvjk = 
5756 K (LAMDA). Therefore assuming that the near-IR 
radiation seen by the HF molecules completely fills the 
sky (/ = 1), the required minimum radiation field tem- 
perature for efficient IR pumping is T > 730 K. This 
corresponds to an energy density vu v ^5x10 erg cm -3 
at 2.5 fj,m, that is, a radiation field >10 6 times that in 
a PDR near to (d = 0.2 pc) an O star (Drainc 2011). 
If / < 1 the radiation temperature, and energy density, 
would be higher. Consequently, if the HF molecules are 
excited by IR pumping in NGC 7130, they should be il- 
luminated by a very intense IR radiation field that can 
be produced only by the AGN. 

6.2.2. Chemical Pumping 

To be efficient, chemical pumping requires that the 
number of HF molecules formed exceeds the number of 
spontaneous decays from the J = 1 to the J = levels 
(see van der Tak 2012). This condition can be expressed 
as 

> AloNlt 



JlHF 



(5) 



where R is the formation rate, Ahf the total HF column 
density, tt-hf the HF density, and N\ the column density 



of the HF molecules in the J=l level. Assuming that HF 
formation and destruction arc balanced: 



R = knYt 2 n F = nn F Qx, 



(6) 



where k = 2.86 x 10 12 cm 3 s 1 is the HF formation rate 
at 50 K, Cd = H7 x 10~ 10 s _1 the HF photodissocia- 
tion rate for xuv? the standard interstellar UV radia- 
tion field of Draine (1978) (Neufeld et al. 2005), x the 
UV radiation field in units of xuv, and tihf and n-F 
are the densities of hydrogen fluoride and atomic fluo- 
rine, respectively. Since tihf + n-p = 2_4f^h 2 > where 
.4f = 2.9 x 10~ 8 is the solar F abundance with respect 
to H (Lodders 2003), Equation 5 can be rewritten as 



CdX 2A F kn H2 
Aw CdX + knn 2 



7Vh 2 > iVi- 



(7) 



For the conditions in a PDR (x = 10 3 , 71h 2 = 
10 4 cm -3 ), the N± calculated before implies a Ah 2 > 
10 24 cm -2 . That is several orders of magnitude 
higher than the H2 column density derived from the 
CO J = 1 — transition (see Section 4.3). Both higher 
gas densities and intenser radiation fields decrease the 
minimum H2 column density. Only a combination of 
dense gas (nn 2 > 10 7 cm~ 3 ) and intense UV radiation 
fields (x > 10 7 ) is compatible with the derived H2 col- 
umn density in NGC 7130. This intense UV radiation 
field suggests that the HF molecules should be close to 
the AGN. 

6.2.3. Collisional Excitation by Electrons 

The rotational levels of molecules are usually excited 
by collisions with H2. However, for molecules with a 
large dipole moment, like HF, the rotational levels can 
be excited by collisions with electrons as well. In gen- 
eral, in the interior of molecular clouds, where most 
of the molecules are found, the electron abundance is 
low (n c /fiH 2 < 10~ 6 ; see Figure 1 of Meijerink et al. 
2006) and the collisions with electrons can be ne- 
glected. However, in some situations, like diffuse clouds 
(Black & van Dishoeck 1991), and molecular clouds sur- 
faces (Wolfire et al. 2010), or XDRs (Meijerink et al. 
2006), the electron density is enhanced. The majority of 
the free electrons in molecular clouds come from ionized 
carbon, with a small contribution from hydrogen ion- 
ized by cosmic rays (Liszt 2012). Therefore, the electron 
abundance can be as high as h c /iih ~ 1.6x10 , that is, 
the gas-phase carbon abundance (Sofia et al. 2004). 

The HF J = 1 — critical density for collisions with 
electrons at 50 K is n c =1.7xl0 4 cm -3 (calculated using 
the LAMDA electron- HF collision strengths). Taking the 
highest n c abundance calculated above, this critical den- 
sity is reached when riH 2 =T0 8 cm~ 3 . This is two orders 
of magnitude lower than the critical density for collisions 
with H 2 , but it is still high. Hence, the HF J = 1 — 
emission line is thermalized only in high-density molec- 
ular clouds. 

Similar to van der Tak et al. (2012), we used RADEX 
to model the HF emission taking into account collisions 
with both molecular hydrogen and electrons. We com- 
puted a grid of models for nn 2 = 10 2 -10 11 cm -3 , T^ n = 
10 — 250 K, and r c /tih 2 = 10~ 4 . We used the electron- 
HF collisional rate coefficients from LAMDA calculated 
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by van der Tak et al. (2012) using the Coulomb-Born ap- 
proximation. The HF J = 1 — line is optically thin for 
a wide range of physical conditions, so the line flux will 
be proportional to -/Vhf in the escape probability approx- 
imation. That is, 

^(ergem" 2 s" 1 sr -1 ) = /C -/Vhf (cm -2 ), (8) 

where K, depends on the density of the collisional part- 
ners, electrons and molecular hydrogen, and on the ki- 
netic temperature. 

For gas densities between ?ih 2 =10 3 -10 4 cm -3 and 
T kin > 50 K, JC ranges between 10 -22 and 10 -21 . The HF 
J = 1-0 flux in NGC 7130 is 1.4xl0- 6 crgcm -2 s -1 sr" 1 
(assuming the SLW beam size, 18", as the source size). 
Accordingly, this corresponds to a beam averaged HF 
column density >10 15 cm~ 2 . Given the F abundance, 
Af = 2.9 x 10~ 8 , the resulting < N H2 > would be un- 
realistically high. For lower Tkin, fC values are lower, 
~10~ 24 -10~ 23 , and the column density would be even 
higher. 

In the high-density limit (riH 2 >10 8 cm -3 ), where 
HF would be thermally excited, JC is ~T0 -18 . There- 
fore, for the observed HF flux in NGC 7130 we obtain 
<7V H F>~10 12 cm- 2 an d <jv H2 >~ 10 20 cm- 2 . It would 

imply ~ 4 x 10 7 M Q of very high density molecular gas 
where C is single ionized. This region would correspond 
to the dark molecular gas (that does not emit in CO), 
which represents about 30% of the total gas in PDRs 
(Wolfire et al. 2010). However typical PDR densities 
are tih 2 ~ 10 4 cm -2 , much lower than the high den- 
sity gas where HF would be excited by electron colli- 
sions. These densities are reached in dense molecular 
cores, but there the UV field is shielded and ionized car- 
bon is not expected. Therefore, the gas where the HF 
J = 1—0 line originates does not seem to be related to the 
star-formation activity in this galaxy. On the contrary, 
X-rays from the AGN can keep a high ionization frac- 
tion over large column densities (Gonzalez- Alfonso et al. 
2013). The importance of X-ray ionization in NGC 7130 
is also supported by the detection of several OH+ tran- 
sitions in its FTS spectrum (see Section 7.2). 

6.2.4. HF Excitation in NGC 7130 

We discussed three possible excitation mechanisms to 
explain the HF ,7=1-0 emission in NGC 7130. These 
mechanisms require high-density molecular gas, in addi- 
tion to strong near-IR or UV emission for the IR and 
chemical pumping, respectively. Although our analysis 
does not favor any of these mechanisms over the rest, 
the physical conditions of the gas are not similar to those 
expected in star- forming regions (PDRs). Consequently, 
the HF J = 1 — emission in NGC 7130 appears to be 
related to the AGN activity. 

7. OTHER MOLECULES 
7.1. Water 

We detected ortho- and para-water (o- and P-H2O) 
transitions in four of our galaxies: two lines in 
UGC 05101, one in NGC 3227, four in NGC 7130, and 
one in NGC 7582 (sec Table 3). All the transition 
are detected in emission. These water lines are weaker 
than the CO lines present in the SPIRE/FTS spectra 



of these Seyfert galaxies. This is also the case in the 
starburst galaxy M 82 and the Seyfert 2 NGC 1068 (see 
Kamenetzky et al. 2012; Spinoglio et al. 2012). On the 
contrary, in the ULIRG Arp 220 water lines are more 
intense than the mid- J CO lines (Rangwala et al. 2011). 

To model the excitation of water using radiation trans- 
fer models we would need to detect more water tran- 
sitions in each galaxy (see e.g., Spinoglio et al. 2012; 
Gonzalez- Alfonso et al. 2012); therefore we did not at- 
tempt this detailed analysis. 

IR pumping is an efficient mechanism to excite water 
molecules. In particular the absorption of far-IR pho- 
tons at 75.4 /im leads o-water molecules to the 32i en- 
ergy level (Gonzalez- Alfonso et al. 2010). In UGC 05101 
and NGC 7582, the o-H 2 3 2 i - 3 i2 line at 1163 GHz is 
the brightest water line in the observed frequency range. 
Thus, in these two galaxies, IR pumping may be im- 
portant for the excitation of water. On the contrary, in 
NGC 3227 and NGC 7130 this line is not detected, hence 
the contribution of IR pumping for water excitation in 
these galaxies could be lower. 

7.2. OH+ 

The formation of OH + is in general associated 
to X-ray /cosmic-ray ionizations, although it is also 
formed in the transition region from ionized to molec- 
ular gas in PDRs (e.g., Sternberg & Dalgarno 1995; 
Hollenbach et al. 2012). This molecule is detected 
in the spectra of galaxies, in absorption in Arp 220 
(Rangwala et al. 2011; Gonzalez- Alfonso et al. 2012) and 
M 82 (Kamenetzky et al. 2012), and in emission in the 
Seyfert galaxies NGC 1068 (Spinoglio et al. 2012) and 
Mrk 231 (van der Werf et al. 2010). Its detection is in- 
terpreted as an indication of XDR excitation in these 
galaxies. 

Three OH + lines are detected in emission in the 
SPIRE/FTS spectrum of NGC 7130 suggesting the pres- 
ence of an XDR in this galaxy. The hyperfine structure 
of these lines is not resolved with the spectral resolution 
of SPIRE/FTS; therefore we measured the sum of the 
fluxes of the hyperfine transitions (Table 4). 

The Einstein coefficients of the hyperfine transitions 
were obtained from the Cologne database for molecu- 
lar spectroscopy (CDMS; MMer et al. 2005). They were 
averaged for each fine structure level: A u i = Ylij Qi-^ijl 
gi, where i and j correspond to each hyperfine upper 
and lower levels respectively. Then, the transition prob- 
abilities for the 909, 971, and 1033 GHz OH + transitions 
are 1.57xl0 -2 , 1.82xl0 -2 , and 2.11xl0 -2 s -1 , respec- 
tively. 

Using Equation 3 we obtained that the beam averaged 
column densities for NGC 7130 are 1.7xlO n , 1.7xl0 n , 
and 1.6xl0 11 cm" 2 for the upper levels of the 909, 971, 
and 1033GHz lines (assuming an 18" beam). 

8. CONCLUSIONS 

We presented the sub-millimeter spectra of eleven ac- 
tive galaxies observed with SPIRE/FTS. Our study is 
focused on the spectral lines tracing the cold and warm 
molecular gas phases in these galaxies. Specifically the 
CO ladder from J up = 4 to 13 lies in this spectral range, 
as well as the two fine structure [C I] lines at 492 and 
809 GHz. For a few galaxies we also detected HF, H 2 0, 
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and OH + lines. In addition, we analyzed new IRAM 
30 m observations of the CO J = 1 — and J = 2 — 1 
transitions for three of these galaxies. The main results 
are summarized as follows. 

1. We modeled the CO SLED (from J up = 1 to 
J up = 7 — 12) of six Seyfert galaxies in our 
sample using the escape probability approxima- 
tion. For five of them the kinetic temperature and 
H2 density of the warm molecular gas are simi- 
lar (n H2 ~10 3 - 2 -10 3 ' 9 cm- 3 and T kin -300-800 K). 
This warm molecular gas seems to be the same 
gas traced by the mid-IR H 2 S(l) 17.03/im rota- 
tional line. The heating mechanism of the warm 
gas is not determined since PDR, XDR, and shocks 
model can be combined to explain the observed CO 
SLEDs. However it is likely related to the star- 
formation activity since the Leo is related to the 
IR luminosity but not to the AGN luminosity in 
these Seyfert galaxies. 

2. We also used a radiative transfer model to inter- 
pret the [C I] emission. The contribution from the 
warm molecular gas to the [C I] emission seems to 
be small in these galaxies. Conversely, this emis- 
sion can arise from cold (Tki n < 30 K) and dense 
(?ih 2 > 10 3 cm~ 2 ) molecular gas where the [C I] 
lines are thermalized. The C° abundance with re- 
spect to H 2 is 3xl0~ 5 -8xl0~ 4 ; therefore the C°/ 
CO abundance ratio is 0.1-2, comparable to that 
observed in star burst galaxies. The [C I] 492 GHz/ 
CO J = 1 — ratio ranges between 15 and 85 in 
our galaxies. It is consistent with that expected 
in PDRs, although for NGC 4388 XDR excitation 
could be possible. 

3. In two galaxies we detected the HF J = 1 — 
transition at 1232 GHz. In UGC 05101 it is de- 
tected in absorption, and implies a column den- 
sity A HF > (1.3 ± 0.3) x 10 14 cm- 2 . In NGC 7130 
the HF transition is observed in emission. We pro- 
posed three excitation mechanisms, near-IR pump- 
ing, chemical pumping, and electron collisions, that 
can produce the HF J = 1 — emission. From our 
analysis the three mechanisms seem to be plausi- 
ble. They require very dense molecular gas, and, 
for the IR or chemical pumping, strong near-IR or 
UV fields. These conditions suggest that the HF 
emission in NGC 7130 is associated with the AGN 
activity. 

4. We are not able to apply any radiative model to 
the water emission. However the detection of the 
o-H 2 3 2 i — 3i 2 transition at 1163 GHz suggests 
that IR pumping (due to strong far-IR continuum) 
could play an important role in the excitation of 
water molecules in some of these galaxies. 

5. In NGC 7130 we detected three OH + lines in emis- 
sion (l o -0i at 909 GHz, l 2 -0i at 972 Ghz, and U 
0i at 1033 GHz). We derived a column density of 
~d.7x 10 11 cm -2 for the upper levels of these three 
transitions. These detections suggest the presence 
of an XDR in NGC 7130. 
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